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■ Summary. We study the dynamics of 2D and 3D barred galaxy analytical models, 

Oh' focusing on the distinction between regular and chaotic orbits with the help of the 

Smaller ALigment Index (SALI), a very powerful tool for this kind of problems. 
5^ . We present briefly the method and we calculate the fraction of chaotic and regular 

' orbits in several cases. In the 2D model, taking initial conditions on a Poincare 

, iHtPy) surface of section, we determine the fraction of regular and chaotic orbits. 

In the 3D model, choosing initial conditions on a cartesian grid in a region of the 
(x,z,py) space, which in coordinate space covers the inner disc, we find how the 
fraction of regular orbits changes as a function of the Jacobi constant. Finally, we 
outline that regions near the {x,y) plane are populated mainly by regular orbits. 
I The same is true for regions that lie either near to the galactic center, or at larger 

^\ ■ relatively distances from it. 

rn 

. 1 Introduction 

OO 

The dynamical evolution of galactic systems depends crucially on their or- 
bital structure and in particular on what fraction of their orbits is regular, or 
, chaotic. Thus to permit further studies, it is essential to be able to distinguish 

;_j ' between these two types of orbits in a manner that is both safe and efhcient. 

This is not trivial and becomes yet more complicated in systems of many de- 
grees of freedom. A summary of the methods that have been developed over 
the years can be found in [1]. 

In the present paper we use a method based on the properties of two 
deviation vectors of an orbit, the "Smaller ALingment Index" (SALI) [2]. It 
has been applied successfully in different dynamical systems [2, 3, 4, 5, 6, 7, 
8, 9, 10, 11, 12], frequently also under the name Alignment Index (AI) [13, 
14, 15, 16, 17] and has been shown to be a fast and easy to compute indicator 
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of the chaotic or ordered nature of orbits. We first recall its definition and 
we then show its effectiveness in distinguishing between ordered and chaotic 
motion by applying it to a barred potential of 2 and 3 degrees of freedom. 
Recently, a generalization of the SALI, the "Generalized ALignment Index" 
(GALI) was introduced by Skokos ct al. (2007) [18], which includes the full 
set of the k initially linearly independent deviation vectors of the system to 
determine if an orbit is chaotic or not. 

2 Definition of the Smaller ALigment Index - SALI 

Let us consider the n- dimensional phase space of a conservative dynamical 
system, which could be a symplectic map or a Hamiltonian flow. We consider 
also an orbit in that space with initial condition P(0) = (xi (0), a;2(0), a;„(0)) 
and two deviation vectors Wi{0), t(?2(0) from the initial point P(0). In order 
to compute the SALI for a given orbit one has to follow the time evolution 
of the orbit itself, as well as two deviation vectors wi{t),W2{t) which initially 
point in two different directions. At every time step the two deviation vectors 
wi{t) and W2{t) are normalized by setting: 

and the SALI is then computed as: 

SALI{t) = min{\\wi{t) + waWH , \\wi{t) - W2{t)\\} ■ (2) 

The properties of the time evolution of the SALI clearly distinguish be- 
tween regular and chaotic motion as follows: In the case of Hamiltonian flows 
or dimensional symplectic maps with n > 2, the SALI fluctuates around a 
non-zero value for regular orbits [2, 3]. In general, two different initial devi- 
ation vectors become tangent to different directions on the torus, producing 
different sequences of vectors, so that SALI does not tend to zero but fluctu- 
ates around positive values. On the other hand, for chaotic orbits SALI tends 
exponentially to zero. Any two initially different deviation vectors tend to 
coincide in the direction deflned by the nearby unstable manifold and hence 
either coincide with each other, or become opposite. 

3 The model 

A 3D rotating model of a barred galaxy can be described by the Hamiltonian 
function: 

H ='^ipl+Pl+ vl) + V{x, y, z) - Qb{xpy - ypx)- (3) 

The bar rotates around its z axis, while the x axis is along its major 
axis and the y axis is along its intermediate axis. The Px,Py and Pz are the 
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Fig. 1. Agreement between the results of the Poincare surface of section (PSS) and 
SALI. a) The Poincare surface of section for the 2D Ferrers model and H = —0.335. 
b) Regions of different values of the SALI for 50,000 initial conditions on the {y,Py)- 
plane integrated up to i = 5000 for the same value of the Hamiltonian. The light grey 
colored areas correspond to regular orbits, while the dark black ones to chaotic. Few 
gray and dark gray colored points correspond to "sticky" orbits. Note the excellent 
agreement between the two methods as far as the gross features are concerned, as 
well as the fact the SALI can easily pick out small regions of instability which the 
PSS has difficulties detecting. 

canonically conjugate momenta. Finally, V is the potential, f2i, represents the 
pattern speed of the bar and H is the total energy of the system in the rotating 
frame of reference (Jacobi constant). The corresponding equations of motion 
are: 



x=px + ^hy, y=Py-QbX, z=p:,, (4) 

dV ^ dV ^ dV 

Px = --g^ + ShPy, Py = --Q^ - ihPx, Pz = ~~Q^- 

The equations of the evolution of the deviation vectors and the calculation 
of the SALI are given by the corresponding variational equations. 
The potential V of our model consists of three components: 

1. A disc, represented by a Miyamoto potential [19]: 

(5) 

^Jx•'+y■' + {A + ^/z' + B■'f 



where Md is the total mass of the disc, A and B are the horizontal and 
vertical scale lengths and G is the gravitational constant. 
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2. A bulge, which is modeled by a Plummer sphere whose potential is: 

GMs 

= =======, (6) 

^/x^ + z-' + ej 

where Cg is the scale length of the bulge and Ms is its total mass. 

3. A triaxial Ferrers bar, the density p{x) of which is: 

P) 

I ,m>l 

where pc = ^if^^^ is the central density, Mb is the total mass of the 
bar and 222 

"^^ = ^ + T9 + ^. a > 6 > 0, (8) 
with a, b and c being the semi-axes. The corresponding potential is: 



where 



A^{u) = {a^+u){b^ + u){c^+u), (11) 

n is a positive integer (with n = 2 for our model) and A is the unique 
positive solution of: 

m\X) = 1, (12) 
outside of bar (m > 1) and A = inside the bar. 

This model has been used extensively for orbital studies [20, 21, 22. 23, 24, 25] 
and we will refere to it hereafter as the Fcrcrrs model. We adopt the following 
values of parameters: G = 1, f2b = 0.054 (54 a = 6, 6 = 1.5, c = 0.6, 

A = 3, B = 1, e, = 0.4, Mb = 0.1, Ms = 0.08,"m^''= 0.82. The units we use, 
are: 1 kpc (length), 1 Myr (time) and 2 x IO^^Mq solar masses (mass). The 
total mass G{Ms + Md + Mb) is set to be equal to 1. 



4 Results in the 2D and 3D Ferrers model 

The 2D Ferrers model is a subcase of the general 3D one and it can be de- 
scribed by the Hamiltonian equation (3) by setting {z,Pz) = (0,0): 



H = -^{pI +pI) + V{x,y) - Qb{xpy - ypa;). 



(13) 
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Fig. 2. a) Percentages of regular (light gray bar), chaotic (black bar) and sticky 
orbits (two intermediate bars - gray and dark gray baxs) for the 2D Ferrers model 
from a mesh of 50,000 orbits on the (y, Py) plane of Fig. 1 and final time of integration 
t — 5000. b) The percentages for 3D, using 50,000 orbits and the same classification. 
The initial conditions are taken as described in the text. 

Fixing H = —0.335, x = 0, we chose 50,000 initial conditions on the {y,py)- 
plane of the Poincare surface of section (PSS), while px = H{x,y,py) and 
we calculated the final values of the SALI for t = 5000. We were able to 
detect very small regions of instability that can not be visualized easily by 
the PSS method. In Fig. 1, on the left panel, we have plotted the PSS for this 
Hamiltonian value. On the right panel we attributed to each grid point a color 
according to the value of the SALI at the end of the evolution. The light grey 
color corresponds to regular orbits and to the areas that host them while the 
black color represents the chaotic ones. The intermediate colors between the 
two extremes represent the so-called "sticky" orbits, i.e. orbits whose nature 
is chaotic but need more time to show their behavior (weak-chaotic orbits). 
The distinction between them is done by measuring SALI a,t t = 5000: orbits 
with SALI< 10~^^ correspond to strongly chaotic, orbits with SALI> 10~^ 
to regular, while orbits with IQ-^ <SALI< lO"* and lO'^^ <SALI< lO'^ 
are "sticky" . In Fig. 2a we present the percentages of orbits from this PSS, 
according to this classification. We find that ~ 37, 3% of the orbits on this 
PSS are regular. 

For the 3D case of the model we used a sample of 50,000 orbits, equally 
spaced on a cartesian grid in the {x,z,py) space, with x G [0.0,7.0], z G 
[0.0,1.5] and Py e [0.0,0.45], while {y,Px,Pz) = (0,0,0). In this way, we 
attempted to create initial conditions that could support the bar. These orbits 
could be parented mainly by the Xi tree, i.e. the Xi family and the XiVi i = 1, ... 
families that bifurcate from it and extend vertically well above the disc region 
[21]. Although these initial conditions cover all the available energy interval, 
they are not spread uniformly over it. Note also that a few of these have an 
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Fig. 3. Percentages of regular orbits as a function of: the Jacobi constant (panel a), 
the absolute value of their z coordinate, <| z |>, meaned over their orbital evolution 
(panel b), their initial spherical radius, Rsphericai (panel c) and the spherical radius, 

< Rsphericai >, meaned over their motion (panel d). 

energy value beyond the escape energy, and we dismiss them. These, however, 
constitute less than 0.5% of the total, so that they influence very little our 
study and the statistics. In Fig. 2b wc present the corresponding percentages 
of regular, chaotic and sticky orbits. We find that, using these initial conditions 
and parameters, the phase space of the model is dominated by chaotic orbits, 
since 77.54% of them are chaotic. 

In order to check the variation of the percentages of regular and chaotic or- 
bits as the energy of the model varies, we did the following: We first sorted the 
energy values for all the initial conditions. Then, wc created 30 energy intervals 
containing equal number of orbits, since our way of giving initial conditions 
does not imply their uniform distribution in the total energy interval. In every 
energy interval wc calculated the percentages of regular and chaotic orbits, 
considering as chaotic all the orbits with SALI< 10~^. In Fig. 3a we show the 
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percentage of regular orbits in each energy interval, as a function of the mean 
energy in that interval. Generally, the percentage of regular orbits decreases 
as the energy increases, but before and after the escape energy (where the Ja- 
cob! constant value is if ~ —0.20) there are two peaks. This non-monotonic 
behavior is related to the appearance or disappearance of stable periodic or- 
bits in the phase space and the size variation of the stability regions around 
them. 

We also attempted to explore the way that regular and chaotic orbits are 
distributed along the ^-direction of the configuration space. Following the 
evolution of each orbit, we calculated the mean of the absolute value of their 
z coordinate (<| z |>). Then, we divided the available <| z |> interval in 
30 slices with equal number of orbits in each one of them. This restriction 
gives us better samples for the estimation of the percentages, implying at the 
same time that these slices are not equally sized necessarily. For every slice 
separately we calculated the fraction of regular orbits and in Fig. 3b we plot 
these percentages as a function of the <| z |> in that slice. It reveals that 
the slices 'near' the (x,y)-plane (<| z \> < 0.35) contain mainly 'regular' 
orbits. Contrarily, slices for larger values of z host mainly chaotic motion. 

Furthermore, we looked at these percentages as a function of the initial 
spherical radius (Rsphericai) and the mean spherical radius (< Rgphericai >, 
meaned over the evolution). Again, dividing in 30 slices the total range of the 
Rsphericai, in & similar manner with the <| z |>, we calculated the percent- 
ages of regular orbits. We plot this for every slice, as a fimction of the mean 
Rsphericai of that slice. We see that the fraction of regular orbits decreases 
strongly with increasing Rsphericai up to Rsphericai < 1.5 where it reaches 
a minimum, while for 1.5 < Rsphericai < 7 this percentage starts increasing 
gradually. This result is in good agreement with the results in Fig. 3d, where 
the horizontal axis corresponds to the value of the spherical radius meaned 
over time during the evolution. 

5 Conclusions 

We used SALI to study the dynamical behavior of Hamiltonian models of 2D 
and 3D barred galaxies. We found that in both cases there is a significant 
amount of chaotic orbits. In the 2D model, we were able to chart a subspace 
of the phase space, to identify rapidly even tiny regions of regular motion and 
measure their percentages. In the 3D model, apart from computing the global 
percentages of regular and chaotic orbits, we calculated these percentages as 
a function of the energy and found that low values of the energy are mainly 
dominated by 'regular' orbital motion. We also followed the distribution of 
the chaotic and regular orbits in the configuration space and we foimd that 
orbits which lie near the (x, 2/)-plane with relatively small mean deviations 
in 2;-direction are generally regular. Finally, we monitored the variation of 
their percentages as a function of their initial spherical radius and their mean 
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spherical radius. We find that the fraction of regular orbits is dominant in 
regions near the center, as well as at relatively larger distances from it. 
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